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1. Introduction 



A recent proposal to obtain d = 2n dimensional lattice chiral fermions from a domain 
wall in d + 1 dimensions [1] has generated a lot of interest. The basic observation is 
that if one introduces a domain wall in an odd dimensional theory of a Dirac fermion by 
introducing a mass parameter which changes sign across a hyperplane of codimension one 
(the "domain wall"), massless states occur that are bound to this domain wall. These states 
therefore can be interpreted as to live in one dimension less, if one takes the size of the mass 
parameter to be of the order of the cutoff (the spread of the bound state wavefunctions is of 
order of the inverse mass in the direction orthogonal to the hyperplane). Moreover, these 
states carry a definite chirality, which depends on the overall sign of the mass parameter. 

Kaplan realized that the fact that the odd dimensional Dirac theory is entirely vector- 
like in nature, without any chiral symmetry appearing at the level of the action, can be 
used to implement this idea on the lattice. Because the action is vectorlike, a Wilson mass 
term can be introduced without breaking any symmetries of the action. This Wilson term 
eliminates the well known species doublers that are generated on the lattice. If we denote 
the coordinates of the d dimensional space by x, and that of the extra d+ 1-st coordinate 
by s, a local {i.e. single site) mass term is added which is positive for s > and negative 
for s < 0. Kaplan showed that in this way the above described phenomenon carries over 
to the lattice, and that a chiral fermion emerges bound to the domain wall defect located 
at s = 0. 

In a sense the theory of a massless chiral fermion constructed in this way is not en- 
tirely free of species doublers. If we consider the theory in a finite but large volume with 
(anti)periodic boundary conditions, an antidomain wall defect has to occur somewhere, 
and a chiral fermion with chirality opposite to that of the massless fermion bound to the 
domain wall occurs at the antidomain wall. In a finite volume the wavefunctions of these 
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two modes have a small overlap (approximately exponentially suppressed as exp{—mL), 
where m is the absolute value of the mass parameter, and L is the distance between the 
domain wall and the antidomain wall), and a tiny Dirac mass couples the two modes, which 
however become massless and decoupled in the infinite volume limit. This phenomenon 
was shown to be quite general {i.e. independent of the precise choice of boundary condi- 
tions) by Narayanan and Neuberger [2]. An alternate, equivalent construction employing 
free boundary conditions on a odd dimensional half space was proposed and discussed by 
one of us [3]. 

This theory was first studied in the presence of background gauge fields in order to 
assess the role of the anomaly. In the original paper [1], it was argued that a Goldstone- 
Wilczek current flowing off the domain wall should exist to reconcile the apparent conflict 
between the anomalous structure of a single even dimensional chiral fermion and the exact 
gauge invariance of the odd dimensional theory, much like the mechanism discussed by 
Callan and Harvey in the continuum case [4]. This was indeed confirmed in a series of 
papers [5,6,7]. In particular, it was shown in ref. [6] that a Goldstone-Wilczek current is 
generated in the region where the combined mass matrix (consisting of the local domain 
wall mass and the Wilson mass terms) is not positive definite, carrying anomalous charge 
between the domain wall and the antidomain wall in a way consistent with the chiral 
zeromode spectrum at those defects [8]. 

The situation becomes more complicated when one considers dynamical gauge fields. 
Obviously, one would like to introduce gauge fields that do not couple both to the ze- 
romodes living at the antidomain wall and also to those at the domain wall, since this 
would render the resulting gauge theory vectorlike. The task is therefore to couple gauge 
fields in such a way that they do not couple to the fermion modes living at the antidomain 
wall. A concrete implementation has been proposed in which the gauge fields are confined 
to a "waveguide" around the domain wall, excluding a region around the antidomain wall 
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[9,10]. Within this waveguide, the gauge fields are taken to be d dimensional, i.e. the gauge 
fields are independent of the direction orthogonal to the waveguide (the s coordinate), and 
all link variables in this extra direction are equal to one [2,9]. (For a different approach, 
where the extent of space in the s-direction is kept strictly infinite, see refs. [11]. See also 
ref. [12].) 

By introducing a waveguide, new defects are necessarily introduced into the theory, 
consisting of the boundaries of the waveguide. These are d dimensional hyperplanes situ- 
ated somewhere "halfway" between the domain wall and the antidomain wall. With the 
specific choice of gauge fields as described above, gauge invariance is explicitly broken at 
these boundaries. This is remedied by introducing a d dimensional scalar field that lives 
only at the two boundaries, restoring gauge invariance in a way similar to the way fermion 
mass terms are made gauge invariant in the electroweak standard model. The Yukawa 
coupling between this scalar field and the fermion field is a free parameter of the model. 
It is important to point out that within the context of the waveguide model, such scalar 
fields are unavoidable. They are essentially the gauge degrees of freedom, which, due to the 
explicit breaking of gauge invariance, couple to the fermions. Moreover, their dynamics is 
necessarily nonperturbative, since their fiuctuations are not controlled by a small param- 
eter (such as the gauge coupling, which only controls the transverse degrees of freedom) 
[13,14]. 

The key question is now of course whether these new defects will lead to new fermionic 
zeromodes. Any such new zeromodes bound to the inside of a waveguide boundary will 
couple to the gauge field, and presumably destroy the chiral nature of the theory. In fact, 
by confining the gauge field to the waveguide, no Goldstone-Wilczek current can reach the 
antidomain wall anymore, and one may legitimately worry whether new chiral zeromodes 
are needed at the waveguide boundaries in order to absorb the anomalous charge carried by 
the Goldstone-Wilczek current between the domain wall and the boundary in the region 
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with nonpositive definite fermion mass matrix. However, it was argued in ref. [10] (which 
from now on we will refer to as i), that this is not necessarily the case, and that the 
existence or nonexistence of new fermionic zeromodes at the waveguide boundaries is a 
dynamical question. 

In order to answer this question one needs to map out the phase diagram and spectrum 
of the theory. This was attempted in /, in which the waveguide model without gauge 
fields, but with the scalar field, was studied analytically and numerically as a function of 
the Yukawa coupling y and the scalar hopping parameter k (the scalar field was taken to 
be radially frozen). It was found that always new zeromodes appear on the waveguide 
boundary in the region of nonpositive definite mass matrix, providing "mirror" modes to 
the chiral modes living on the domain wall in a one-to-one correspondence. 

However, as described in /, the numerical simulations were notoriously difficult for 
large values of the Yukawa coupling, due to a bad signal to noise ratio for the fermion 
propagators in that region, and it was hard to arrive at any definite conclusion for large 
y. In the present paper we return to the question of the fermionic zeromode spectrum 
at the waveguide boundaries for large Yukawa coupling. It is well known that in general 
new phases occur for large values of the Yukawa coupling of lattice scalar-fermion models. 
Those phases can be distinguished from those at small values of the Yukawa coupling by the 
fermionic zeromode spectrum. In particular, it is in principle possible that no zeromodes 
occur at all at large y. This would be an important ingredient for constructing a genuinely 
chiral lattice gauge theory from domain wall fermions. 

Here, we develop an expansion suitable for investigating the large y region of the phase 
diagram. Again, we set the gauge field link variables equal to one, but keep the scalar 
field at the waveguide boundary. This allows us to find the zeromode spectrum for large 
y. We do find a new symmetric phase which is not connected to the symmetric phase at 
small values of y. However, the zeromode spectrum in fact turns out to be very rich, and 
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quite similar to that of the Smit-Swift [13,15,16] model for loea/i; Yukawa coupling: all the 
doublers of the zeromode at the domain wall reappear at the waveguide boundaries! 

In section 2 we define the waveguide model restricting ourselves to the case d = 4, 
and recast it in a way amenable to an expansion in 1/y. In the next section we find the 
zeromode spectrum for y = oo, and in section 4 we estimate the location of the phase 
transition between the symmetric and broken phases in the K-y plane for large y. We then 
briefly consider the zeromode spectrum at large but flnite y. In the last section we present 
our conclusions. 

2. The model 

The waveguide implementation of the domain wall fermion model is defined by the 
action 

^domainwall = ^^^^ " ^(^) + + ^ - W + rUs)^ s 

seWG s0A^G 

+ [*si"i?*s+l + *s+l^L*s]- 5^^5*5 (1) 

In this equation we have suppressed the four dimensional coordinates, and only indicated 
the fifth coordinate explicitly, is the four dimensional symmetric nearest neighbor 
difference operator, and w denotes the four dimensional Wilson mass term. D^{U) and 
W{U) are gauge covariant versions of these, rus is the domain wall mass, which is taken to 
be positive for s > and negative for s < 0. (We can take m(0) = for instance.) V is a 
radially frozen scalar field which takes its values in the gauge group G. The waveguide WG 
consists of all lattice points with s'^ < s < sq, and the domain wall at s = is inside the 
waveguide {i.e. Sq < and sq > 0). With periodic boundary conditions, there will be an 
antidomain wall outside the waveguide, where rus again changes sign. We have chosen the 
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Wilson parameter r = —1, and the projection operators Pr[l) ^Ii{L) ~ + (~)75)- 



For full details, see /. (Here we have chosen some conventions different from those in /, 
where for instance we chose r = 1.) This model contains a righthanded massless fermion 
bound to the domain wall and a lefthanded massless fermion bound to the antidomain 
wall, irrespective of the value of y (c/. i). 

In the present paper we are interested in the fermionic zeromode spectrum at the 
waveguide boundaries, and we will restrict ourselves to a simpler model in which only the 
waveguide boundary between s = sq and s = sq + 1 is present, with no other defects. We 
imagine all other defects to be very far away. Also, as in /, we will choose U = 1, i.e. we 
will turn off the gauge fields. Including a hopping term for the scalar field, the action for 
this situation is 



S = Sp + Sb, 



(2) 



with 




X,Sy^S0 





X 



and 




(4) 



where 




(5) 



□ 



xy — 



'x-n,y 




(6) 



We will always choose < m < 1. 
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This action is invariant under the transformations 

*a;,s 9^x,s, s < sq, 

■^x,s^h-^x,s, s>so + l, (7) 
V^gVh\ 

where g and h are elements of G. 

For y small and in the symmetric phase {i.e. k < Kc) one massless lefthanded fermion 
exists just inside the boundary, and a righthanded one just outside, with wavefunctions 
which fall off exponentially away from the boundary. In the broken phase, these two chiral 
modes pair up into a massive Dirac fermion with mass proportional to yv with v the 
vacuum expectation value of the field V . Combined with the zeromodes at the domain 
wall and the antidomain wall, and turning on the gauge fields, this model results in a 
vectorlike gauge theory for small y {cf. T). 

We would now like to investigate the situation for large y. An expansion in 1/y can 
be developed by a slight rewriting of the action, eq. (3). First, we will relabel the fermion 
variables 

Xs = *s, s< so, 
XRso = *i?so' XLso = ^LsQ+h (8) 
Xs = *s+l, s > so, 

where r^l — Pr,L^ ^ ^^c. Note that -0 is a four dimensional fermion field, i.e. independent 
of s. Rescaling 

V' ^ ^V', (9) 

the fermionic action becomes 

SF^Y.^^sHs + Xs(^{s){u+m-l)xs) (10) 
s 
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+ {XLsXRs+l + XRs+lXLs) + '^L^'^R + '^rVHl 

s 

+ Va(V'(n + m - l)Xso + Xso(° + ™ - + 

with 

a = - (11) 

y 

and 

a(s) = l-5,,,o. (12) 
In eq. (10) we again suppressed all x dependence. 

3. The fermion spectrum on the boundary for y ^ oo 

For y ^ oo, i/j decouples, and the Dirac equation for Xs{p) — ^~^^^Xx,s is 

{p)Xs{p) + a(s)(m - 1 - F{p))xs{p) + PrXs^i{p) + ^LXs-l(p) = 0, (13) 

with ^ (p) = ^n^n^vapn and F{p) = X^/xCl ~ cosp^). We wish to look for solutions of 
this equation (in Minkowski space) which also satisfy the four dimensional massless Dirac 
equation, and which have a definite chirality. For such solutions 

XRs+l{p) = a(s)(l - m + F{p))xRs{p), (14) 
XLs-liP) = - TO + F{p))xLs{p)- 

These equations have normalizable solutions bound to the waveguide boundary due to the 
fact that a(so) = 0. The first equation has solutions inside the waveguide (i.e. XRs — foi" 
s > sq) for F{p) > m. Likewise, the equation for XL has solutions outside {i.e. XLs = 

for s < Sq) for F{p) > m. 

This leads to a rich spectrum of chiral zeromodes at the waveguide boundary. Let 

TTAeSuO, (15) 
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£ = {(0, 0, 0, 0), (tt, tt, 0, 0), ... , (tt, tt, tt, tt)} 

= {(tT, 0, 0, 0), (0, TT, 0, 0), ... , (0, TT, TT, Tt)}. 

The four dimensional Dirac equation has relativistic (continuum) solutions for p = tt^ + p 
withp — > 0. The chirality gets flipped if tt^ G O, i.e. xr{p) is lefthanded whenp ^ tt^ e O 
[17]. The condition F{p) > m excludes the modes around p = 0, and results in a spectrum 
of 7 righthanded and 8 lefthanded zeromodes inside the waveguide boundary, and the 
mirror reflection of this outside the boundary. With the single righthanded zeromode on 
the domain wall and the lefthanded one on the antidomain wall the zeromode spectrum is 
nonchiral both inside and outside the waveguide. The fact that the spectrum outside is a 
mirror copy of the spectrum inside is in accordance with a symmetry of the action: 

Xx,s 74XPa;,2so-s' 
ipx lAi^Px^ (16) 

where Px = {—xi, —X2, —x^, X4). 

The other boundary of the waveguide can be treated in exactly the same way, with 
now however m replaced by — m, since this boundary is located on the other side of the 
domain wall. The condition for the existence of zeromodes becomes F{p) > —m on both 
sides of this boundary, which is satisfied for all p. It follows that there are 8 righthanded 
and 8 lefthanded zeromodes both inside and outside this boundary for y ^ 00. 

We will later need the x-fermion propagator for y — > 00, which can be found as the 
inverse of the Dirac operator 

{[i^ (p) + a(s)(m - 1 - Fip))]Ss,s" + PrSs+1,s" + PlSs-1,s") = Ss,s'- (17) 

We find, for s' = sq 

(0) _ -z/^(p) + ^(p)-6(p) 

^sMp)- sHp) + mHp) °' 
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where 



= 1 - m + > 0, Vp, (19) 

z{p) = >1, Vp, 

M^{p) = 1 - b{p)/z{p) > 0, Vj9. 

Let us check that this propagator gives rise to the same zeromode spectrum as we found 
from the Dirac equation. For p — tt^ we have b{p) — 2n + 1 — m, with n equal to the 
number of components of tt^ that are equal to tt. This leads to 2; = 1/(1 — m) for n = 
and z = 2n + 1 — m for n ^ {1,2, 3, 4}, and hence the propagator has a pole at p = tt^ 
when TTj[ ^ 0, but not for p = 0. 

For the other waveguide boundary b{p) — l+m + F{p), and in this case z = 2n+l + m 
for all n, so that the propagator has massless poles at all tt^, in accordance with what we 
found earlier. 



4. The phase diagram for large y 



In the limit ?/ — > 00 the fields V' ^ind V decouple from each other, and the phase 
diagram of the model is dictated by S]^, in particular by the bosonic hopping parameter 
K. For G = U{1) or SU{2) there will be a second order phase transition at k = Kc- (The 
fermionic part of the action depends on V even for y cxd, but the fermion determinant 
is independent of V in this limit.) The theory consists of a scalar field in the broken or 
symmetric phase, and free fermionic modes, described by the field x- 

We would now like to see how this changes if we turn on the coupling a. The effective 
action for V can be obtained by integrating out the fields and x ^in expansion in a. 
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Let us first integrate out Xi which can be done exactly, resulting in an efi^ective action for 



with 



where 



^i^=Yl {^LxVxi^Rx + i^RxVxi^Lx) (20) 

X 

+ a ^ (i^Jxyi^y - bP{n +m- l)]xS^^^^^^^^y^^^^ [(□ + m - l)^/']^) , 
xy 

4"!..).(...) = //*'^""'*»°?..W, (21) 



'\ (22) 



(27^) 

To order the effective action for V obtained by integrating out all fermion fields then 
becomes 

SeffiV)^SB{V) (23) 
- logdet [5^,2/ + c^iviPu + VyPi) {l?x,y - [(□ + m - l)^S,lo(n + ™ - l)]x,j/) 

-2^-gt,vj..)E/e-(^-^)(--»)v(.)v(.) ,.j^i,.(^),.,;;t,.(^) 



+ 0(a^) 



xy fj, 
4^ 



(recall s^(p) = sinp^). 

There will be a similar contribution from the waveguide boundary on the other side of 
the domain wall, with the only difference that m be replaced everywhere by — m. 

This effective action contains long range self-interactions for the field V, due to the 
fact that the integrand of the last term has massless poles at p = tt^ 7^ (and for all tt^ 
for the terms with m —>■ —m). It is well known that a meanfield approach gives the correct 
qualitative phase diagram for ck = 0. While this is in general not true for more complicated 
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interactions, we expect it to hold for a. small enough. If we apply meanfield, we obtain an 
estimate for the location of the critical line in the k-ck plane (valid for small a) : 

1 



K + -a'[I{m) + I{-m)]+0{a'') = Kc, (24) 
with 

f s\p)h\p) , 1 f s\p)h\p) 

JpsHp)+M\p)^ 2jp[s^{p)+M\p)r ^ ^ 

which is a function of m through h and M. As an example, the value of / is 187.4(3) 

for m = 0.5 and 371.3(6) for m = —0.5. The large values of I{m) reflect the enhanced 

coupling of the fermions to the V field (c/. the factors h{p)). 



5. Fermion spectrum for large but finite y 

We can integrate out the field '4> in an expansion in a = 1/j/ in order to obtain an 
effective action for x and V: 

Sx=Y. ^^s^Xs + Xs(^{s){u + m - l)xs) + (XsPrXs+I + Xs+iPlXs) (26) 
s s 

oo 

- aY^Xsoi^ + rn - I) J](-a)"[(Ftp^ + FPi)^]"(Ftp^ + vPl){u + m - l)x.o 

n=0 

s,s' X 

+ «^ E[^«o(° + ™ - ^)]x{v1Pr + V:cPL)^xy{viPR + FyPi)[(n + m - l)xso]y 

x,y 

To order ct, replacing F by its expectation value f 1, the Dirac equation reads (in momen- 
tum space) 

{p)Xs + a(s)(m - 1 - F{p))xs + J^iJXs+l + PlXs-1 -av{m-l- F{p)f6s,soXso = 0. (27) 
Effectively, a{s) of eq. (12) no longer vanishes at s = sq. Instead, 



a(so) = av{l-m + F{p)), (28) 
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and there are no zeromode solutions for v ^ ^. In the broken phase, the zeromodes pick 
up a mass of order v/y. 

We see that in the broken phase, we effectively obtain a nonzero value for a(so) (eq. 
(28)). This is only possible in the broken phase. The symmetry cq. (7) acts on Xso ^ 

XRsq ^ 9XRso'> XLsq ^ ^X.Lsqi 

and this symmetry has to be broken in order to have a nonzero value for a(so). Thus, 
a righthanded mode inside the waveguide can mix with a lefthanded mode outside the 
waveguide only if the symmetry eq. (7) is violated, i.e. only in the broken phase. More- 
over, no righthanded and lefthanded modes which are both inside the waveguide can mix, 
because no lefthanded mode inside the waveguide carries the same lattice momentum as 
any righthanded mode inside the waveguide, and lattice momentum is exactly conserved. 
(See the last section for further discussion of this point). 

We will illustrate this by showing that the massless poles which exist at y — > oo (c/. 
eq. (18)) remain massless to second order in a = 1/y, and that therefore the massless 
spectrum of the theory remains the same for large but finite y. For G = U{1), to order 
CK^, the inverse propagator in the symmetric phase becomes 

^s,Ap) = -io^y {p)6s,soSso,s'b\p)iL + K{p)), (29) 

with 

L =< VxV^+^ > (30) 

and 

/^('-'''^fSo-^WifW. (31) 

D(j)) is the scalar propagator. The L term comes from a diagram with one order vertex 
(c/. eq. (26)), and the K{p) term comes from a diagram with two order a vertices. Since 
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vanishes for all p = tt^ (eq. (15)), the massless pole structure of the propagator is not 
changed by the o? correction. (Note that the function K goes to a constant for v ^ 
in the symmetric phase.) We believe this result to be valid to all orders in a. 

Let us conclude this section by considering the effective action for eq. (20). This 
effective action is exact, and therefore should describe all zeromodes at the waveguide 
boundary for all values of the Yukawa coupling j/, as long as V' couples to the field Xi 
i.e. for CK 0. The term linear in a contains the y — > 00 propagator for the field x, 
and therefore the zeromodes at the waveguide boundary at large y are present in the 
effective action S^. For small y, it was shown in / that there is only a zeromode at the 
waveguide boundary for tt^ = 0, and this should also follow from S^. This is indeed the 
case: by rescaling the field t/^ in a way appropriate for small values of y (i.e. undoing the 
transformation eq. (9)), and choosing y = 0, we see that the ijj propagator for small y is 
defined by the order a term in eq. (20) to lowest order in y. In momentum space, this 
propagator reads 

This expression only has a pole at p = 0, and none at p = tt^ with tt^ 7^ (c/. end 
of section 3), and describes one massless Dirac fermion at the waveguide boundary, in 
accordance with the results of / at small y. 

6. Conclusion 

In this paper we continued the investigation of the waveguide implementation of the 
domain wall fermion model [9,10]. In / we argued that the question as to whether this 
model can yield a chiral gauge theory is a dynamical question which cannot be decided on 
the basis of simple anomaly arguments. Therefore, a complete investigation of the phase 
diagram is needed in order to see what the massless fermion spectrum of the model is, as 
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a function of the parameters of the model. However, in / this analysis was left incomplete, 
because of difficulties with the numerical exploration of the phase diagram for large values 
of the Yukawa coupling. 

Here we showed that the action for this model can be reformulated in such a way that 
a systematic strong coupling expansion in the Yukawa coupling can be performed. This 
allowed us to examine the strong Yukawa coupling region of the phase diagram analytically. 

In the limiting case a = (j/ — > oo) the fermions decouple from the scalar field, and 
we obtained a free fermion theory which can be solved exactly. We found a very rich 
zeromode spectrum. While the zeromode spectrum at the domain wall and antidomain 
wall is unchanged, in a sense all their doubler modes show up localized at the waveguide 
boundaries. The complete massless fermion spectrum inside as well as outside the waveg- 
uide is vectorlike. The gauge fields, which live only inside the waveguide, will couple to all 
the massless fermions inside, and the result is a vectorlike gauge theory. The situation is 
in fact worse than at weak Yukawa coupling, where inside the waveguide only one massless 
mirror mode emerged due to the interaction of the fermions with the scalar field. 

The phase diagram consists of symmetric phases at weak and strong Yukawa coupling 
(for small enough values of the scalar hopping parameter), separated by a phase in which 
the gauge symmetry is spontaneously broken. We verified by an explicit calculation that 
the massless spectrum remains unchanged to first order in in the symmetric strong 
coupling phase. While we did not provide a rigorous proof, there can be little doubt that 
this result extends to all orders in a^. There is an interesting relation between the present 
model and the Smit-Swift model [18]. The massless fermion spectrum at strong Yukawa 
coupling is similar to that of the Smit-Swift model at weak Yukawa coupling, and vice 
versa. 

The question arises whether a more general framework could lead to qualitatively 
different results. For example, by a suitable modification of the scalar action, it might 
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be possible to arrange for the existence of an antiferromagnetic condensate which is g- 
and /i-invariant. In that case, lattice momentum would be conserved only modulo tt (and 
not modulo 27r) and mixing between righthanded and lefthanded modes from different 
corners of the Brillouin zone would be possible inside the waveguide boundary. Without 
fine tuning, we would expect that only a single lefthanded mode would survive inside the 
boundary. But since lefthanded and righthanded modes would pair to form massive states, 
the full massless spectrum inside the waveguide would remain vectorlike! 

There is a good reason to believe that the phenomenon described in this example is 
completely general. The point is that we start (at a = 0) with a well-defined theory that 
has a relativistic, albeit vectorlike, low energy spectrum. The question is whether one can 
achieve a chiral spectrum by making a continuous change of parameters which leaves the [g 
and h) symmetries unbroken. Assuming the model maintains a consistent continuum limit 
throughout this process (which should be the case if we add only finite range operators to 
the lattice action), we expect that the low energy spectrum can be correctly described by 
a corresponding change of parameters in a relativistic effective lagrangian. But then the 
only way to decouple massless fermions is to pair them into massive Dirac fermions. Thus, 
since we started with a vectorlike spectrum, we will end up with a vectorlike spectrum. 

We conclude that, although the phase diagram is richer than suspected in ref. /, the 
waveguide implementation of the domain wall fermion model does not lead to a lattice 
chiral gauge theory. 
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